Conductivity, shot noise, and hot phonons in bilayer graphene 
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We have studied electrical conductivity a and shot noise of bilayer graphene sheets at high bias 
voltage Vds- As a function of bias, we find linear increase of a which is leveled off above Vds ~ 0.2 
V. We assign the saturation of a to the creation of optical/zone boundary phonons. Excess noise 
Fano factor F is found to increase first with bias and reach a maximum at Vds — 0.1 V, above 
which F decreases as V d ~~" with 7 = 0.25 at the charge neutrality point. Using mean-free-path 
type of analysis we find a good agreement between the calculated and measured optical phonon 
temperatures. 

PACS numbers: PACS numbers: 73.63.-b, 73.23.-b, 73.50.Fq, 73.50.Td 



Monolayer and bilayer graphene are highly promising 
candidates for new semiconductor materials [l[ . In mono- 
layer form, it is gapless as its cone-like conduction and va- 
lence bands touch at two inequivalent Dirac points, where 
the density of states vanishes. The conductivity at the 
Dirac point, however, remains finite [lLyJ and the shot 
noise emulates that of a diffusive wire Bilayer 
graphene has attracted large interest partly due to the 
fact that an energy gap could be opened by chemical dop- 
ing [?J or by applying external electric field In gen- 
eral, the basic transport properties of bilayer graphene 
at the charge neutrality (or "Dirac") point (CNP) are 
very close to monolayer graphene: their shot noise is ex- 
pected to be equivalent [9( while the conductivity scales 
with the number of layers. When tuned away from the 
CNP, more involved behavior is predicted for the bilayer 
Q , though experimental verification of the differences is 
still missing, presumably due to disorder [lfj| . 

At large voltage bias, ballistic transport in both mono- 
layer and bilayer graphene is expected to break down 
due to emission of optical phonons (OP), as has been ob- 
served in single- walled carbon nanotubes (SWNT) [ill ]. 
In monolayer graphene, the hig h bias behavior has been 
investigated by Meric et al. [12j who observed a tendency 
to current saturation due to OP emission. Electron- 
phonon (e-ph) scattering is detrimental for the speed and 
for the transconductance of graphene devices, whereas for 
the shot noise OP emission may contribute in a beneficial 
way by suppressing current fluctuations. Here we concen- 
trate on bilayer graphene and report the first studies of 
shot noise in it. We show that at large bias voltages the 
shot noise is, indeed, decreased due to phonon emission. 

Our theoretical treatment of resistance in a graphene 
bilayer is based on simple mean-free-path arguments, 
as has been done for example in Ref. [TJ for SWNTs. 
The resistance R, dependent on drain-source voltage 
bias Vds, is defined through the current I(Vd s ) = 



We assume R to have the form [14| R{V C 



ds 



sample, R ba ii{Vds) = R q /M e ff(V ds ) is the conductance 
for a ballistic junction, with R q — h/2e 2 the quantum re- 
sistance and M e ff(Vd s ) an effective number of transverse 
modes within the transport window, and R c is an addi- 
tional "contact resistance". Scattering is taken into ac- 
count with an effective momentum-relaxation mean free 
path l e ff(Vd s ). Using Matthiessen's rule l e ff(Vd s ) = 

,(Vds)) j where the individ- 
ual mean free paths are due to impurities, OP emission, 
and OP absorption, respectively; acoustic phonon scat- 
tering is neglected due to its weakness [15|, [l6| . Here k mp 
also depends on the carrier density n (or gate voltage V g ) 



J'imp ^op,em(^ s ) ^op,abs\ 



as well as the impurity density 



M eff {V ds ) too 



depends on V g via the Fermi energy. 

We assume that all the relevant OPs have the same 
energy Kl w 0.2 eV, and that they all have the 
same occupation number N op = N(T op ). Here T op is 
the corresponding temperature, and we define N(T) — 
l/(exp(HQ/k B T) - 1). The bath temperature T = T Q 
determines the equilibrium occupation No — N(Tq). For 
emission, we assume l op ,em(Vds) = L ]3t~\ + le ~°P N° + +i ' 
where the first term is the distance needed for the elec- 
trons to reach the threshold energy Ml for OP emission 
and the second term is the mean free path after the 
threshold has been exceeded 11, For absorption, we 
take l p,abs(Vds) = le-op N ^ +1 ■ It is expected that when 
\V ds \ -> 00, then N op -> 00, and so l op . em , l op ,abs -> 0. 

To find out the stationary-state Vd s -dependence of N op 
or, equivalently, of T op , a heat balance equation must 
be solved self-consistently with the current. We employ 



Rball{Vds)—J-f 



R c . Here L is the length of the 



similar heat flow schematics as for SWNTs 13. 
bottleneck in the heat transfer from the electrons to the 
substrate is assumed to be due to the decay of OPs 
into acoustic or substrate phonons (see also Ref. Il9f ). 
and the latter are assumed to remain at T = To. The 
power that heats the OP system is P = Pj ou i e — Pdiff- 
Here Pj ou ie — I 2 {R — R c ) is the Joule power and 
Pdi ff = 2n 2 {k B /e) 2 {T 2 - T 2 )/{R - R c ) is the power 
transferred from the electrons to the leads. Here T„ is 



2 



an average electron temperature. In stationary state 



P — P 

N c nn(N(T e ) 

N q )/t op - 



N op )/t, 
Here P P _ r 



where we approximate P e -op 
_ op and P op -ac = N c htt(N op 



and 



(-^op— ac £md T op ^ ac ) 



are an electron-OP (OP-acoustic phonon) power trans- 
fer and relaxation time, respectively. Assuming that 



-o P 1~op—ac 



we find T„ 



T op , where T op is given 



by N op = N + P/(N c hfl/T op _ ac ). The parameter N c 
is an effective number of occupied modes on the various 
OP branches. A value N c < LW/A C , where W is the 
width of the sample and A c the unit cell area, models 
the nonuniformity of the OP distribution over the phase 
space jilj. We note that some of the above approxima- 
tions may only be justified for Vds > 0.2 V, where both 
e-op and electron-electron (e-e) scattering are strong. 

The lower two electronic energy bands of a bilayer are 
given by E(k) = ±{^Jt 2 _ + (2hv F k) 2 -t_i_)/2, where 1 1 = 
0.4 eV is the interlayer coupling and vp = 1.1 TO 6 m/s [1|. 
This yields the number of ballistic propagating transverse 
modes at energy E: M{E) = ^W-^E 2 + t±\E\. At 
Vds t± /e also the higher bilayer bands should begin to 
contribute, but we neglect them, assuming the scattering 
from phonons to be very strong there. Using this we es- 
timate M eff {V ds ) = ^ J^+I e ^»l M(E)dE. The effec- 
tive Fermi energy fj, > is used to model the presence of 
charge puddles and evanescent modes, such that a resid- 
ual carrier density n* remains also at the CNP (2^|. We 

2 / dn 



(Cg]V g 



V g ,CN P \leYI%, 



use the form [i = (n*) 2 
where C g is the specific gate capacitance and = ^i v -i 
is the density of states. 

Shot noise has turned out to be very useful in the stud- 
ies of graphene @, @j. At the zero frequency limit, shot 
noise is given by the correlation function of current fluc- 
tuations 5I(t): Si = J dt(SI(t)5I(0)) [23]. Typically, the 
strength of shot noise is given by the Fano factor, defined 
as 5" = Si/2e(I). For bilayer graphene, theoretical cal- 
culations give either | @] or 1 — |- i.e. very close to 
what has been observed at the CNP in monolayers 0, Q • 

At high bias, inelastic OP scattering reduces 5/ in gen- 
eral. The role of inelastic scattering on the shot noise in 
diffusive conductors has been investigated by Nagaev [25 1 
and by Naveh et al. [26[ . They find a power law decrease 
of the Fano factor # ~ f3~ 2/5 where j3 = L/l ln > 10, 

kn = (lap,em( V ds) + Cplafcs (Vds))' 1 ■ Thc result is based 

on a solution of the Boltzmann equation and the expo- 
nent of the power law reflects the e-ph collision kernel 
which is sensitive to the dimensionality of the sample, to 
the nature of the phonons, as well as to the details of 
the thermal transport. For ID phonons, the Boltzmann 

— 1/3 

equation yields g" oc V ds ' . 

The above result can be obtained using a heating 
model [27] which assumes that one can define a local 
temperature T{x) along the sample, independent of the 
transverse dimension of the sheet. Moreover, this model 




FIG. 1: (Color online) I(Vd s ) curve measured between the 
terminals H3 and H4 on a 0.3 x 1.3 fim 2 sheet at 4.2 K at the 
CNP. Optical interferogram of the sample is illustrated in the 
right inset. To make the sheet more visible, digital contrast 
enhancement has been applied. The upper inset displays our 
typical bilayer Raman spectrum with four lines fitted. 



yields a simple relation for the Fano factor 

2k B T e 



3 



eV ds 



(1) 



where T e = J Q L T(x)dx is the average electronic tem- 
perature. This result is valid when either e-ph or e-e col- 
lision rate is large and local equilibrium is approached. In 
graphene both e-op and e-e scattering are strong and bias 
dependent (lq : the e-e mean free path is ~ L already at 

v ds - 0.1 V. 

Inelastic scattering by e-e interaction can also lead to 
an increase of the Fano factor. This takes place if the heat 
carried away via the electron-phonon coupling is small 
compared with the (electronic) heat diffusion along the 
graphene sheet to the metallic contacts. With uniform 
resistivity and Wiedemann-Franz law, the above integral 
yields T e = ^V ds , i.e. # = ^ [23] for hot electrons in 
local thermal equilibrium. Now, if heat leaks out from 
the system via e-ph coupling, the ratio ksTe/eVds goes 
down and g decreases. Hence, we expect to have two 
opposite tendencies for $(Vd s ), but at very large bias, 
the e-ph coupling will dominate and reduce the noise. 

We have measured three bilayer graphene samples hav- 
ing L = 300 nm and widths W = 1.0, 1.3, and 1.6 /Ltm, 
respectively. The samples were contacted using Ti/Al/Ti 
sandwich structures with thicknesses 10 nm / h j 5 nm 
where h was varied over 50 — 70 nm (10 nm of Ti is 
the contact layer). A set of 0.3-0.6 ^m wide contacts, 
2 — 6 per sheet, were patterned using e-beam lithog- 
raphy. The strongly doped Si substrate was used as a 
back-gate, separated by 250 nm of Si02 from the sam- 
ple. All the data were measured at To = 4.2 K using a 
two-lead configuration. Differential AC-conductivity was 
recorded at frequency / = 32 Hz while the excess noise 
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FIG. 2: (Color online) Differential conductivity a for bilayer 
sample with L = 0.3 /im and W = 1.3 /im. versus Vds- The 
gate voltage values are indicated in the image. The CNP is 
located at V g = V g ,cNP = —38 V (lowest curve). The circles 
show the dependence of a(Vds = 0) on V g scaled to the Vds 
axis, as discussed in the text. 
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FIG. 3: (Color online) a) Result of our model for the CNP 
(lower curve, V g = -38 V in Fig.[2} and for V g = V 9 ,cnp +38 
V using parameters given in the text, b) Electron temperature 
T e as deduced from measured shot noise using Eq. {TJ (o), and 
OP temperature T op = fcBln ((jv ftn |_i)/jv — 7 fr° m our model with 
the same parameters as in (a) (solid line). 



F = (Si (I) — Si(0))/(2eT) was integrated over the range 
/ = 600 - 900 MHz. The relation between F and the 
true Fano factor J is given by 



eVds 
2k B T 



coth 



eVds 
2k B T 



1 / 



eVds 
2k B T 



(2) 



L dl 

W dV d3 

We measured 



For further techniques, we refer to Ref. H|. 

The I{Vds) curve of our 0.3 x 1.3 ^m 2 sample is illus- 
trated in Fig. [TJ It displays non-linear behavior [28| at 
small drain-source voltage Vds] the nonlinearity depends 
on the gate voltage V g and it is largest near the CNP. The 
non-linear behavior can be observed more clearly in Fig 
[5] which displays the differential conductivity a 
of the same sample at a few values of V g 
the minimum zero-bias conductivities ctq = 6.5, 5.5, and 
4.7e 2 //i for our samples of size 0.3 x 1.3, 0.3 x 1.6, and 
0.3 x 1.0 yum 2 , respectively. The enhancement from the 
universal value 2 x 4e 2 /irh w 2.5e 2 /h suggests the 
presence of impurities in these samples. 

A nearly linear increase in a is seen at Vds — 40 — 200 
mV in Fig. [2j the slope of the increase depends weakly 
on the distance from the CNP which is located at V g = 
Vg,CNP = — 38 V. The increase in a is an indication that 
M e ff grows as the bias window is increased. Below 40 
mV, the conductance displays interference phenomena 
(not shown), reminiscent of the Fabry-Perot oscillations 
observed by Miao et al [2j . The growth in a slows down 
above 0.2 V and a saturation is reached above Vds = 
0.5 V. At larger Vds, even a down turn of a is found. 
This indicates a tendency towards saturation of current, 
although it's clearly much weaker than typically seen for 
carbon nanotubes llj or recently in monolayer graphcne 
flil I29I . However, our samples are much shorter than 
HQ, 



those of Refs. 1 



and the saturation is known to be 



weak for short nanotubes (30j . 

Figure [3^, illustrates conductivity results of our model 



calculation using the parameters: Ml — 0.18 eV, C g — 
100 aF//im 2 (see below), R c = 70 Vt per contact, 
imp = 9 • 10 10 cm~ 2 , l e ^ op = 15 nm, T op _ oc = 5 ps, 
N c = 0.2(LW/A C ), and n* = 9 ■ 10 11 cm" 2 . Here the 
V g dependence of cr(Vd s — 0) fixes rii mv and n*, while 
SWNT values are used for r op _ ac [2l| and l e - op 
also Ref. [13|). Note that N c and r op _ ac only appear in 



[l| (see 



the combination N c MI/t 



It should also be noted 



that the independence of the contact resistance on V g was 
recently shown by Russo et al. [3l| to be a good approxi- 
mation. Furthermore, the fit is done bearing in mind the 
OP temperature deduced from shot noise measurements 
(see below). The agreement between the model and the 
experimental results is good. 

The finding, implicit in the above fitting, that a(Vd s ) 
is related to the effective number of modes in graphcne 
can be given further support by comparing the measured 
Vds- and Vg-dependencies of a. We find linear scaling 

a (Vds =0,V g = ^v)=a (V ds = V, V g = 0). Fitting 

this to the data (see Fig. [2]) yields the value C g = 100 
aF//zm 2 used above, which is close to the theoretical 
parallel-plate value C g = 140 aF//xm 2 . 

The measured excess noise Fano factor F is illustrated 
in Fig. [4] for our 0.3 x 1.3 /zm 2 sample on log-log scale. 
Owing to its definition, F increases first linearly in the 
cross-over regime between thermal and shot noise. The 
dashed curve, fitted to the low-bias data at the CNP 
(V g = —38 V) indicates the behavior obtained from Eq. 
([2]) using $ — 0.28. One can see that at intermediate 
bias, Vds = 50 — 100 mV, the Fano factor grows beyond 
the low-bias value, indicating overheating of electrons in 
the sample. However, the observed value is well below 
$ = V3/4 so that we do not reach pure hot electron 
regime. Using small bias fits Vds < 40 mV, we find # = 
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tion and a slight decrease of conductivity due to phonon 
scattering; T op given by our model agrees well with the 
values deduced from shot noise measurements. 

We thank A. Bachtold, Y. Hancock, A. Harju, T. 
Heikkila, A. Karkkainen, M. Laakso, P. Pasanen, S. 
Russo, E. Sonin, and P. Virtanen for fruitful discus- 
sions. This work was supported by the Academy of Fin- 
land, EU contract FP6-IST-021285-2 (CARDEQ), the 
NANOSYSTEMS project with Nokia Research Center, 
and the CIMO exchange grant KM-07-4656 (MW). 



FIG. 4: (Color online) Excess noise Fano factor F — (Si (I) — 
Si(0))/(2el) vs. bias voltage Vds at three values of V g indi- 
cated in the figure (from top to bottom). The dashed curve 
is calculated using Eq. <(2j) for 5 = 0.28. 



0.23 and $ = 0.2 for our other two samples at the CNP. 
These values fall somewhat short from the values of # = 
1 — 2/ 71", 1/3 calculated theoretically for bilayer graphene 
9, 24]. However, disorder may decrease shot noise and 
improve the agreement (see e.g. Refs. l32l [33h . 

All the traces in Fig. [4] display a maximum near Vds = 
0.08 V, nearly independent of the gate voltage value al- 
though the actual maximum value for F decreases from 
0.31 at the CNP down to 0.23 far away from it. We note 
that the position of the maximum coincides with the in- 
flection point of a (Vds) (see Fig. [2j. Above 0.1 V, all the 
curves display a power law decrease F cx V^" 7 with an 
exponent 7 = 0.25 — 0.13. The decrease in F already 
at Vds ~ 0.1 V suggests that there are active OP modes 
around 100 meV as well. We have also checked our 0.3 x 
1.3 /1m 2 sample with a 15 nm thick hafnium oxide ALD 
coating on it and obtained different exponents, which fits 
the expectations as the Coulomb screening and the OPs 
become both modified. 

Taken together, all evidence indicates short inelastic 
mean free path at Vds ^ 0.1 V. Consequently, we argue 
that local T(x) is well defined at large bias and that shot 
noise can be employed to determine the electronic tem- 
perature T e according to Eq. |T]). Fig. [3b displays T e at 
the CNP obtained from the data in Fig. 01 For compar- 
ison, we have plotted the OP temperature T op from our 
model. The agreement is good, and we conclude that 
both the electrons and the OPs heat up to 1000 K in our 
measurements. 

We have presented data on high bias electric transport 
on bilayer graphene sheets with disorder. At Vds > 0.1 
V, we observe a decrease of Fano factor in a power- 
law fashion, which can be interpreted as shot noise re- 
duction by optical/zone boundary phonons with energy 
Ml ~ 0.2 eV. Using a model for resistance, based on sim- 
ple mean- free-path type of arguments, we find an optical 
phonon mean free path which is close to that determined 
in SWNTs. At bias voltages of 0.5 - 1 V, we find a satura- 
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